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Abstract 

If primordial scalar and tensor perturbation spectra can be inferred from 
observations of the cosmic background radiation and large-scale structure, then 
one might hope to reconstruct a unique single-field inflaton potential capable 
of generating the observed spectra. In this paper we examine conditions under 
which such a potential can be reliably reconstructed. For it to be possible 
at all, the spectra must be well fit by a Taylor series expansion. A complete 
reconstruction requires a statistically-significant tensor mode to be measured 
in the microwave background. We find that the observational uncertainties 
dominate the theoretical error from use of the slow-roll approximation, and 
conclude that the reconstruction procedure will never insidiously lead to an 
irrelevant potential. 



PACS numbers: 98.70.Vc, 98.80.Cq 



1 Inflation and Perturbations 



Inflation produces metric perturbations, which are presently the most plausible cause 
of the observed temperature fluctuations in the cosmic background radiation (CBR) 
and which may act as the seeds for structure formation. During inflation the Hubble 
radius H^ 1 (H = a/a is the expansion rate) must have increased more slowly than the 
scale factor a, while during the radiation and matter eras it increased as a 2 and a 3//2 , 
respectively. The dynamics of the evolution of the Hubble radius (or, equivalently, the 
evolution of the expansion rate H) during inflation is usually modelled by assuming 
that the energy density is dominated by scalar-field potential energy. The number of 
dynamical degrees of freedom associated with the evolution of the scalar-field potential 
energy (and hence the Hubble radius) is unknown. If there is only one relevant degree 
of freedom, that degree of freedom is called the inflaton, and the scalar-field potential 
is called the inflaton potential. 

Inflation models where a single inflaton field slowly evolves under the influence of 
the inflaton potential are called single-field slow-roll models. Such models are very 
attractive because of their simplicity, because they arise naturally in a host of particle 
physics models, and because other more complicated models can often be expressed 
in terms of an effective single-field slow-roll model. 

The assumption that only a single scalar field is dynamically relevant can, at least 
in principle, be tested by a series of consistency relations which represent relations 
between the scalar and tensor perturbations produced in single-field models [[I], 0. 
Inflation models where more than one scalar field is dynamically important are phe- 
nomenologically much more complicated, and almost certainly must be studied on a 
case-by-case basis. In this paper we shall restrict our discussion to the single-field 
case, where the scalar potential V((f>) is the only piece of information to be specified. 

Reconstruction of the inflaton potential (see Ref. || for a review) refers to the 
process of using observational data, especially microwave background anisotropies, to 
determine the inflaton potential capable of generating the perturbation spectra in- 
ferred from observations [Q . Of course there is no way to prove that the reconstructed 
inflaton potential was the agent responsible for generating the perturbations. What 
can be hoped for is that one can determine a unique (within observational errors) 
inflaton potential capable of producing the observed perturbation spectra. The re- 
constructed inflaton potential may well be the first concrete piece of information to 
be obtained about physics at scales close to the Planck scale. 

As is well known, inflation produces both scalar and tensor perturbations, which 
each generate microwave anisotropies (see Ref. U for a review). As reviewed in the 
next section, if V(<p) is given the perturbation spectra can be computed exactly in 
linear perturbation theory through integration of the relevant mode equations. If 
the scalar field is rolling sufficiently slowly, the solutions to the mode equations may 
be approximated using the slow-roll expansion [|, ||, |[]. The standard reconstruction 
program makes use of the slow-roll expansion, taking advantage of a calculation of the 
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perturbation spectra by Stewart and Lyth [|7|], which gives the next-order correction 
to the usual lowest-order slow-roll results. 

The biggest hurdle for successful reconstruction is that many inflation models pre- 
dict a tensor perturbation amplitude well below the expected threshold for detection, 
and as we will see, even above detection threshold the errors can be considerable. If 
the tensor modes cannot be identified a unique reconstruction is impossible, as the 
scalar perturbations are governed not only by V(<f>), but by the first derivative of V(<f>) 
as well. Knowledge of only the scalar perturbations leaves an undetermined integra- 
tion constant in the non-linear system of reconstruction equations. Another problem 
is that simple potentials usually lead to a nearly exact power-law scalar spectrum 
with the spectral index close to unity. In such a scenario only a very limited amount 
of information could be obtained about high energy physics from astrophysical ob- 
servations. 

Recently, Wang et al. || suggested an alternative danger to reconstruction. They 
studied the possibility that the potential might be sufficiently complicated that the 
slow-roll expansion breaks down, and a full integration of the mode functions is nec- 
essary to calculate the perturbation spectra to the accuracy expected by the next 
round of satellite observations by map and Planck. 

It is, of course, clear that any perturbative procedure such as the slow-roll expan- 
sion has some limit to its validity. But here, the breakdown of the slow-roll expansion 
is particularly troublesome; firstly because of the technical difficulty that there is no 
known method of computing yet higher-order corrections, and secondly because there 
is a general possibility that for sufficiently complicated potentials the perturbation 
series may not converge at all. However, the key concern for reconstruction is whether 
such a situation will lead to a misconstruction of the inflaton potential. 

In this paper, we shall argue that for the purpose of reconstruction the breakdown 
of the slow-roll approximation is in many senses desirable. First of all, we show 
that there is no practical danger that one can be misled into thinking perturbative 
reconstruction is working when in fact it is not. Secondly, the very failure of the 
perturbative approach implies that extra parameters are required in order to explain 
the observed data set, and hence there is the prospect of obtaining a greater amount 
of information about the inflationary model. Finally, we will demonstrate that the 
theoretical errors from the use of slow-roll equations are sub-dominant to the expected 
observational errors from a satellite experiment such as Planck. 

Wang et al. || described a particular, admittedly rather ad hoc, scalar field poten- 
tial, for which they demonstrated that the result of exact mode equation integration 
differed significantly, by the standards of say the Planck satellite, from the slow-roll 
expansion. We shall largely concentrate on that potential as a means of concrete 
illustration of the very general points we shall make. 
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2 The slow-roll expansion 



Two crucial ingredients for reconstruction are the primordial scalar and tensor per- 
turbation spectra As(k) and Ar(/c) which we define as in Ref. 0. We define slow-roll 
parameters e, rj, and £ as 
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where here and throughout the paper the prime superscript implies d/d<p. As long as 
these parameters are small compared to unity, As(k) and Ar(k) are given by [|7J 
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where H((f>), e(0), and ^(0) are to be determined at the value of when k = aH 
during inflation, and where C = —2 + In 2 + 7 ~ —0.73 is a numerical constant, 7 
being the Euler constant. These equations are the basis of the reconstruction process 
and provide an analytic alternative to integration of the full mode equations. 

The expressions for A 2 s {k) and A?n(k) are in terms of H(<f>) and the slow-roll 
parameters e(0) and rj(4>), but what is usually known is the potential V(4>). Since H 
can be written exactly in terms of V and e: 
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the problem reduces to finding e and rj in terms of V and its derivatives. 

There are two methods to determine H and the slow-roll parameters. The first is 
an analytic approach studied by Kolb and Vadas and by Liddle et al. 0. (The 
results are identical; here we follow the notation of the later paper.) The approach 
involves an infinite hierarchy of potential slow-roll parameters: 
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and so on, where = d n V/dcf) n . 
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The slow-roll parameters of Eqs. fl2.1| - [2.3|) involve an infinite number of potential 
slow-roll parameters, and to obtain an analytic expression it is necessary to truncate 
the expression by neglecting potential slow-roll parameters above a given order. For 
most potentials the higher derivatives become small rapidly, and the expressions can 
be terminated after only a few terms. For instance, to fourth order in H 2 H] 
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and similar expressions can be found in Refs. || [| for e(0), r]((f)), etc. 

In practice, one would calculate the first few potential slow-roll parameters, and 
if the ones involving higher derivatives become small, the series may be truncated. 
However if the higher-derivative slow-roll parameters do not become small, then it is 
necessary to perform a numerical calculation. 

The slow-roll parameters may be calculated numerically quite easily by expressing 
Eq. (|2.6|) as an equation for e: 
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taking the derivative with respect to 0, using the expressions for H' and H in terms 
of e [Eq. (|2.1|) and Eq. (|2.6|)1 to yield a first-order differential equation satisfied by e, 



3). 



(2.10) 



It is easy to show that r\ = e — mpi e'/ ^/I6ne, so using Eq. (|2.10| ) for e' and integrating 
it to find e, rj is found. The slow-roll parameter £ can be expressed in terms of e, 
T], and 7/ as £ 2 



erj — mpir}'Je/47r. For this expression the additional function e" is 
required, which can be expressed as 
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For the calculation of A${k) and Ax{k) there are three levels of ease/accuracy 
one may employ as indicated in Table 1. The simplest, but least accurate, method 
is an analytic approach, using Eq. (|2.8j ) and corresponding ones for H, e, and rj, and 
the slow-roll expressions, Eqs. (^) and ( |2.5| ), for A 2 s (k) and Aj,(k). This will be 
accurate if the potential slow-roll parameters e v , rj v , etc. become small, so that a 
truncation can be performed in the expressions (the analogues of Eq. ( |2.8| )) for the 
slow-roll parameters e and 77, and provided the slow-roll parameters are smaller than 
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Table 1: Different approaches to computing the power spectrum from a potential 

v{4>). 



unity. If e and r] are small enough, the square brackets in Eqs. (|2.4|) and ( |2.5| ) may 
be replaced by unity and the expressions simplify further. 

The most accurate and reliable, but most involved, is a numerical calculation in 
which the slow-roll parameters are calculated numerically and used in a numerical 
integration of the mode equations. 

For some potentials, it may be desirable to perform something we call a semi- 
analytic calculation, which involves calculating e and rj numerically, and using the 
numerical values in the slow- roll expression for A 2 s {k) and A^(k). This will be useful 
if the higher-order potential slow-roll parameters become large but the actual values 
of e, t], and £ are small. Only one numerical integration is required to compute the 
spectrum, whereas the mode equation integration must be done for each k. 

3 The potential 

Wang et al. || consider an inflaton potential of the form 

K(0)=A 4 (l--tan- 1 ^) , (3.1) 
\ it m Pl J 

which is shown in Fig. [I]. We shall take A = 1; results could be rescaled if desired to 
match the COBE normalization. With this potential, inflation is everlasting. Wang et 
al. |H envisage extra physics, such as a hybrid inflation mechanism |TUfl , intervening to 
end inflation at some point, and assume that the field's location at the time observable 
perturbations were formed can be placed anywhere on the potential. As a specific 
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Figure 1: The potential of Eq. |37T . 



example, they choose <ft — — 0.3mpi to be the value of the field when the present 
Hubble radius equalled the Hubble radius during inflation. We adopt their choice for 
the end of inflation. 

Fig. ^] shows the potential slow-roll parameters of Eq. ( |2.7| ), r v and ( v being the 
next two in the hierarchy. Since potential slow-roll parameters become large, the 
truncation necessary to calculate the slow-roll parameters in terms of the potential 
slow-roll parameters is not reliable in the vicinity of the origin, though they have 
chosen the potential carefully so that the violation is not sufficient to even temporarily 
end inflation by driving e larger than unity. 

Given the approximation of linear perturbation theory, the perturbation spectrum 
can be computed exactly taking the inflationary dynamics into account by using the 
mode equation formalism of Mukhanov (Tl|. Numerical solution of this equation 
was described by Grivell and Liddle |12| and Wang et al. The scalar and tensor 



perturbations generated by this potential are shown in Fig. [| We also show the result 
of the best available semi-analytic technique as described above, which uses numerical 
solution of the background field equations as input to the Stewart-Lyth formula and 
which gives an error of more than ten percent at some k. We see that the matter power 
spectrum is very far from the scale-invariant form, featuring a sharp downturn and 
then bending back up again. In fact, this deviation from scale-invariance is already 
too large to be permitted by observational data; for example, in a COBE-normalized 
CDM-like cosmology the perturbations at 8ft. -1 Mpc are only o~g ~ 0.3. It still serves 
as a useful example however. 
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Figure 2: The first six slow-roll parameters. 



Any set of observations will be able to probe only a range of scales, and we shall 
use the capabilities of the Planck satellite as a benchmark. Planck can probe the 
microwave anisotropy multipoles from the quadrupole I = 2 up to around i = 2000. 
The corresponding k range sampled is from k = aoHo to k = 1000aoi/o- Using 
a H = /i/3000Mpc~\ the range of k probed by Planck is shown in Fig. |3|. We 
see that in fact Planck would be able to see more or less to the bottom of the dip, 
but not the subsequent rise. So the part of the spectrum to be probed does not 
contain any very strong features. Fig. |] shows the predicted Ct curves, both exact 
and approximate, calculated using a modified version of the cmbfast code [T3|. The 
slight differences from the equivalent figure in Ref. |§ are presumably from slightly 
different choices of parameters such as the baryon density. The basic slow-roll method 
does extremely poorly. 

To illustrate how to recognize when reconstruction is impossible, we shall later 
consider the case where the scales left the horizon a little further down the scalar 
field potential, corresponding to shifting the spectrum leftwards by a factor of 30 and 
bringing the main feature of the spectrum into the observable range. 



4 The Reconstruction Procedure 



In this section we review the general reconstruction procedure, and extend results for 
high derivatives of V to higher-order in slow-roll than previously given. Perturbative 
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Figure 3: The perturbation spectra obtained by mode equation integration (solid 
curves) compared to the semi-analytic results (dotted curves) and the lowest-order 
analytic results (dashed curves). Since the amplitude of the power spectra are pro- 
portional to the unknown parameter A in the potential, A 2 S and A\ may be adjusted 
up and down by a common factor. Since inflation is terminated by hand, the spectra 
can also be shifted left or right. Shown in the figure are the spectra analyzed in 
subsection 5.1. The spectra analyzed in subsection 5.2 are shifted to the left by a 
factor of 30 so that the dip in A 2 S is near k = 0.02/iMpc -1 . 

reconstruction requires that one fits an expansion, usually a Taylor series of the form 
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to the observed spectrum in order to extract the coefficients, where stars indicate the 
value at k*. The scale fc* is most wisely chosen to be at the (logarithmic) center of 
the data, so we take A;* = 0.01 h Mpc -1 . 

How far the series should be taken is governed by how many of the coefficients can 
be obtained with error bars inconsistent with zero. Planck can in principle measure 
the Ct at the percent level of much of its range, though degeneracies with other 
parameters will weaken the determination of the spectrum somewhat. This sets a 
ballpark figure for how far the series should go, and for almost all known inflation 
models only the amplitude and spectral index are required. 



8 




Figure 4: Predicted microwave anisotropies from the power spectra computed under 
different assumptions summarized in Table 1: from top to bottom they are numeri- 
cal, semi- analytic, and lowest-order analytic, corresponding to the scalar and tensor 
spectra shown in Fig. || 



In general one should similarly fit the tensor amplitude. The relative amplitude 
of tensors and scalars is important for reconstruction. If more information about 
the tensors can be determined than just the amplitude, then the extra information is 
degenerate with information in the scalars and the consistency relations can be tested, 



though it seems extremely unlikely that observations will permit this p], [14], [L5 . 

From these, we reconstruct the potential, using the equations set down in Ref. 
plus similar equations for V" and V"" } This requires the exact relation connecting 
changes in with changes in k 0, 
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The reconstruction equations, using the notation n* to indicate d % n/d In k % \^ and 
R* = A?^(^k^ / A'g (^*) j are 

75 

m^V{^) ~ — A 2 S {K)R*{1 + [0.21R*)} 

m^V'fa) ~ -^^A 2 S (K)R^ 2 {1 - [0.85i^ - 0.53(1 - n*)] } 



1 For the potential we reconstruct the V"" term is negligible. 
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We can assign a bookkeeping parameter to keep track of the order of the terms 
appearing in Eq. (|4.3j) . /2, and 1 — n* is of order 1, and n* is of order j + 1. One can 
see that Eq. ( |4.3| ) is arranged so that the higher-order corrections appear in square 
brackets. Since the fit we use for the scalar power spectrum only goes up to n* , we 



have no information about n* 
expression for V"". 
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hence for consistency we only use the lowest-order 
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. . . are to be determined from observations. Fortunately, parameter 



estimation from the microwave background has been explored in some detail p4], |T5" 



We shall use error estimates for Planck assuming polarized detectors are available, 
following the analysis of Zaldarriaga et al. |I5| . Most analyses have assumed that 
/2, and n* are the only parameters needed to describe the spectra. In a related 
publication [I6"|, we have generalized their treatment to allow the power spectrum to 
deviate from scale-invariance. Including extra parameters leads to a deterioration in 
the determination of all the parameters, as it introduces extra parameter degeneracies. 
Fortunately, for most parameters the uncertainty is not much increased by including 
the first few derivatives of n [|16|] , but we have found that the parameter n itself has a 
greatly increased error bar. If a power-law is assumed it can be determined to around 
An ~ 0.004 1 15]) |16 1 , but including scale dependence increases this error bar by a factor 
of ten or more. Notice that unless one assumes a perfect power-law behaviour, this 
increase in uncertainty is applicable even if the deviation from power-law behaviour 
cannot be detected within the uncertainty. 

From Ref. ||16|| , an estimate of the relevant uncertainties is 

A(/2,)~ 0.004 ; An ~ 0.14 

A{dn/d\nk) ~ 0.04 ; A[d 2 n/d(\n kf] ~ 0.005 . (4.4) 

Note that r ts in Ref. |16[] is 14 times /2,. 
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5 Two Reconstruction Examples 



In this section we consider two reconstruction attempts of regions of the potential in 
Fig. 1. One is successful and the other unsuccessful. We demonstrate that this is 
to be expected and is perfectly consistent with the reconstruction procedure we have 
adopted, indicating that it is virtually impossible to misconstruct the potential. 



5.1 Successful reconstruction 

In this subsection we reconstruct the potential shown in Fig. |l] that led to the spectra 
of Fig. ^. We find that the observational errors dominate over the errors introduced 
by using the slow-roll approximations for the spectra. We do not agree with claims in 
Ref. || that the breakdown of the slow-roll approximation for this potential illustrates 
that reconstruction is not valid for general potentials. We also find that, as expected, 
the shape of the potential is better determined than the amplitude. 

To estimate the errors introduced by the slow-roll formalism, we treat the exact 
power spectra of Fig. |3] as the 'observations', and imagine that we do not know the 
potential from which they derive. We shall later use estimated uncertainties assuming 
the spectra are measured via the microwave anisotropies, but for now we shall imagine 
that the power spectra have been measured exactly in order to assess the errors arising 
from uncertainties in the theoretical power spectrum determination. 

We have fit a Taylor expansion of the above form to the exact spectrum, using k* 
as above, and found the following results: 



lnA 2 s {K) = 3.84 ; #, = 0.014 
dn 



n* = 0.579 



dink 
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-0.052 . (5.1; 



From the estimated observational uncertainties of Eq. ( }4.4| ), we see that all these 
coefficients are successfully determined at high significance and a simple "chi-by-eye" 
demonstrates that the spectrum reconstructed from these data is an adequate fit to 
the observed spectrum. This stresses the point that the more unusual a potential 
is, the more information one is likely to be able to extract about it, though the 
uncertainties on the individual pieces of information may be greater. 

We reconstruct the potential in a region about 0* = —0.22 (the reconstruction 
program does not determine 0*), the width of the region given via Eq. ( |4.2| ). The 
results for the reconstructed potential without observational errors are shown in Fig. 
[5], where it can immediately be seen that the reconstruction has been very successful 
in reproducing the main features of the potential while perturbations on interesting 
scales are being developed. 

Thus far we have assumed perfect observations, but of course in practice the 
measurement of the power spectrum is subject to observational errors. To illustrate 
visually the effect of the observational errors, we use a Monte Carlo procedure where 
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Figure 5: The original potential is compared with the fragment reconstructed from 
Eq. (|4.3|) using the parameter determination in Eq. (|5.1| ). The lower figure is the 
reconstructed m 3 Pl V' /V 3 ^ 2 . The difference between the reconstructed fragments and 
the true potential is due to the use of the slow-roll expansion in reconstruction. 
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Figure 6: A set of ten reconstructed potential fragments compared with the true 
potential. The spread in the reconstructed potentials compared to the reconstructed 
potential in Fig. |^ illustrated our point that the uncertainty in the parameter de- 
termination will overwhelm the error introduced by using the slow-roll formalism in 
reconstruction. 
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we draw 'observed' parameter values from gaussian distributions about the true value 
with the appropriate variance. These observed parameters, which include fully the 
effects of cosmic variance and the modelled instrument noise, are then fed into the 
reconstruction process. Fig. | shows ten such reconstructions about the true value. 
The uncertainty is indicated by the wide spread in the results. Note the width of the 
reconstructed region is correlated with the height of the potential; if the observed -R* 
happens by chance to be low the potential is lower and flatter, so the field rolls more 
slowly as the perturbations are generated. 

The uncertainty is dominated by that of A\{k*); although the gravitational waves 
are detectable in this model, it is only about a three-sigma detection and the error 
bar is thus large. Since the overall magnitude of the potential is proportional to A?n, 
the visual impression is of a large uncertainty. 

Fortunately, information in combinations of the higher derivatives is more accu- 
rately determined. Fig. ^ also shows the reconstruction of V /V 3 ^ 2 with observational 
errors; this combination is chosen as it is independent of the tensors to lowest order. 
Not only is it reconstructed well at the central point, but both the gradient and cur- 
vature are well fit too, confirming useful information has been obtained about not just 
V" but V" as well, which is only possible because of the extra information contained 
in the scale-dependence of the power spectrum. So despite the poor visual impression 
given by Fig. ^ rather accurate information is being obtained about the potential. 

To compare these observational errors with the purely theoretical errors of the 
slow-roll approximation, one can compare to Fig. |5|. From this we conclude that the 
error introduced by using the bessel function approximation in reconstruction is a 
small part of the error budget. 

5.2 An example where reconstruction is impossible 

What we have seen in the example of Subsection 5.1 is that despite slow-roll not 
working very well, and hence generating significant errors in the normal estimation 
of the power spectrum, the reconstructed potential still looks very good, especially 
in the light of expected observational error bars. The main reason for this is that 
although the spectrum is quite far from scale-invariance, it could still be reasonably 
fit by just a few terms in the Taylor expansion. 

The truly awkward situation would be a spectrum which could not be adequately 
fit that way, for then perturbative reconstruction would break down. This is no 
surprise — there is always a limit beyond which a perturbative process cannot be 
taken. An example, alluded to above, would be if inflation ended just a little later 
on the potential we have been considering, so that the spectrum is shifted to the left 
bringing the minimum into the observable range. 

We have examined this case and found that indeed the spectrum cannot be fit 
by a Taylor series to any reasonable order. That reconstruction fails is completely 
obvious, because we can't even start. There is no danger of trying to carry out 
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a reconstruction, and obtaining an answer which has no connection to reality. If 
the observations take this unfortunate turn, then one has a number of options to 
try. Perhaps single-field inflation is not correct at all, and some other theory such 
as topological defects or isocurvature fits better. Or perhaps we do believe that a 
complicated inflaton potential is at work, in which case a non-perturbative technique 
(e.g. fitting the power in different wavebands^) would be more appropriate. Or a 
multi-field inflation model might be responsible, and one could test them on a model- 
by-model basis to see if they were compatible with the data. 

6 Conclusions 

We have found that a complicated potential, such as the example used by Wang et al. 
||, would be a boon to reconstruction of the inflaton potential, as it would provide 
extra information accessible to observations, in the form of scale-dependence of the 
density perturbations. Such potentials are also more likely to have tensor modes at 
a detectable level, which is required for a complete reconstruction to be performed. 

We have analyzed the impact of observational errors on the reconstruction, and 
conclude that any theoretical errors from use of the slow-roll equations are likely to 
be sub-dominant. One might even be able to test this be solving the mode equation 
for the reconstructed potential, though we have not tried to pursue that route here. 

Finally, it has always been clear that there must be some limits to the applicability 
of perturbative reconstruction, although it does a good job even for the Wang et al. 
potential. However, such a failure should be immediately evident from the data, as 
one would not be able to successfully match observations with a truncated Taylor 
series expansion of the power spectra. It would appear, therefore, that there is no 
danger that perturbative reconstruction might appear to be working, but in fact be 
producing an irrelevant potential. 
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